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Abstract 

This  paper  explores  the  possibility  of  improving  the  power  density  of  a  solid  oxide  fuel  cell  (SOFC).  A  three-dimensional  computational  model 
(CFD-ACE  package),  with  the  relevant  sub-models  was  used  for  the  study.  The  performance  of  the  SOFC  was  examined  with  a  thin  wall,  which 
splits  the  inlet  section  and  runs  up  to  half  the  length  of  the  flow  channels.  The  results  obtained  with  this  (thin- walled)  geometry  were  consistently 
better  than  those  obtained  with  plain  geometry  (without  the  thin  wall).  The  polarization  characteristics  of  the  thin-walled  geometry  indicated  that 
the  maximum  power  density  obtained  was  1.18  Wcm~2  at  an  efficiency  of  around  60%.  The  corresponding  values  of  maximum  power  density 
and  the  efficiency  at  which  it  was  obtained  for  a  plain  geometry  were  0.88  W  cm~2  and  50%,  respectively.  The  enhanced  performance  of  the 
thin-walled  geometry  was  attributed  to  a  better  distribution  of  the  reactants  along  the  length  of  the  SOFC.  Studies  were  also  conducted  to  verify 
the  performance  of  the  thin-walled  geometry  over  a  wide  range  of  inlet  mass  flow  rates.  They  revealed  a  superior  performance  of  the  thin-walled 
geometry  compared  to  the  plain  geometry.  At  lower  inlet  mass  flow  rates,  the  difference  between  the  two  in  performance  was  small,  but  at  higher 
inlet  mass  flow  rates  the  difference  in  performance  was  significant. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Following  the  discovery  by  Nernst  in  1899  that  zirconia  can 
act  as  an  oxygen  carrier,  a  practical  solid  oxide  fuel  cell  (SOFC) 
was  developed  by  Baur  and  Pries  in  1937.  A  good  historical 
perspective  of  the  development  of  SOFC  can  be  found  in  [1]. 
Recently,  major  consortiums  like  Siemens-Westinghouse  [2,3], 
BMW  Delphi  [4],  Rolls  Royce  [5]  have  all  identified  the  low 
power  density  of  SOFC  as  an  area  for  improvement,  if  the  future 
SOFC  systems  are  to  become  financially  viable  following  the 
target  set  by  the  Solid  State  Energy  Conversion  Alliance  ( SEC  A) 
[6]. 

A  brief  review  of  some  of  the  past  efforts  to  enhance  the  power 
density  of  a  SOFC  is  presented  here.  Li  et  al.  [7]  have  proposed 
a  new  design  for  the  gas  distributor  in  fuel  cells,  which  improves 
the  power  density  of  a  SOFC.  In  their  design,  a  large  number  of 
small  current  collectors  surrounded  by  a  small  reactive  area  have 
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been  used  to  obtain  the  enhanced  power  density.  They  claim  that 
through  their  new  design  the  access  area  of  reactants  to  the  elec¬ 
trode/electrolyte  layer  has  led  to  an  enhanced  mass  transfer  and 
higher  power  density.  They  have  demonstrated  the  viability  of 
their  new  design  through  experiments  on  PEM  (proton  exchange 
membrane)  fuel  cells,  wherein  they  claim  to  have  obtained  a  40% 
improvement  in  power  density  for  a  single  fuel  cell. 

Sung  et  al.  [8]  and  Takashi  et  al.  [9]  have  reported  obtain¬ 
ing  improved  power  density  in  SOFCs  by  using  a  new  mate¬ 
rial.  Sung  et  al.  [8]  have  modified  the  Ni/yttria-stabilized  zir¬ 
conia  (YSZ)  cermet  anode  by  coating  with  a  samaria-doped 
ceria  (SDC,  Smo.2Ceo.8O2)  sol  within  the  pores  of  the  anode. 
They  report  a  50%  increase  in  power  density  along  with  a 
smaller  degradation  in  performance  over  time.  Takashi  et  al. 
[9]  report  an  increased  maximum  power  density  (from  161  to 
213mWcm-2)  using  an  yttria-stabilized  zirconia  solid  elec¬ 
trolyte  with  a  25wt.%  Ceo.8Gdo.20i.9  (GDC)-containing  Ni 
anode  and  a  15wt.%  Mn02-containing  Lao.8Sro.2Mn03  cath¬ 
ode. 

Bharadwaj  et  al.  [2]  and  Lu  et  al.  [10]  have  performed 
computational  studies  on  the  flattened  tubular  SOFC  designed 
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Fig.  1 .  The  current  density  variation  along  the  axial  length  of  the  SOFC  at  y  =  0 
and  z  =  0  at  a  voltage  of  0.73  V. 


by  Siemens-Westinghouse  Corporation.  In  these  studies,  it  is 
reported  that  a  flattened  SOFC  delivers  a  higher  power  den¬ 
sity  than  the  conventional  tubular  SOFC.  Kim  et  al.  [11]  have 
developed  an  anode-supported  flat-tube  SOFC  by  an  extrusion 
process  and  they  report  of  having  obtained  225  mW  cm-2  at 
750  °C.  Similar  power  densities  have  also  been  reported  by 
Shao  et  al.  [12]  while  working  on  a  thermally  self-sustaining 
single-chamber  micro  SOFC  operating  on  propane,  oxygen  and 
helium  mixture  (helium  being  the  inert).  A  detailed  account  of 
studies  that  discuss  the  fabrication  of  alternative  geometries  for 
SOFCs  to  improve  the  power  density  can  be  found  in  Singhal 
and  Kendall  [13], 

In  this  paper,  a  different  approach  is  taken  to  explore  the  pos¬ 
sibility  of  enhancing  the  power  density  of  a  SOFC.  Fig.  1  shows 
the  current  density  variation  along  the  length  of  a  co-flow  SOFC 
obtained  from  computations  (detailed  explanations  of  how  this 
was  obtained  are  presented  later  in  this  paper)  at  a  voltage  of 
0.73  V.  The  fuel  and  the  air  enter  the  SOFC  system  at  x  =  0  and 
leave  the  system  at  x  =  40  mm.  The  current  density  as  seen  from 
Fig.  1,  decreases  from  the  inlet  section  to  the  exit.  Similar  trends 
have  been  observed  by  Camapanari  and  Iora  [14]  while  carrying 
out  simulations  on  planar  SOFC.  Experimental  results  on  the 
current  density  measurements  performed  by  Matti  et  al.  [15]  on 
a  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  also  indi¬ 


cate  a  similar  variation  along  the  axis  of  the  fuel  cell.  This  aspect 
can  also  be  understood  from  the  basic  Butler- Volmer  equations 
(refer  [16]).  These  equations  suggest  that  the  EMF  developed 
is  strongly  dependent  on  the  partial  pressures  of  the  reactants 
and  products.  It  decreases  with  the  decrease  in  the  partial  pres¬ 
sure  of  the  reactants.  As  one  moves  from  the  inlet  section  of  the 
SOFC  towards  the  outlet  in  a  co-flow  arrangement,  the  reactant 
mass  fraction  or  its  partial  pressure  continuously  decreases  due 
to  the  consumption  of  the  reactants  and  production  of  products 
in  the  reaction.  The  interplay  of  the  above-mentioned  factors  is 
depicted  in  the  paper  by  Yakabe  et  al.  [17].  The  decrease  in  the 
partial  pressure  of  the  reactants  due  to  the  consumption  of  the 
reactants  and  production  of  products  in  the  reaction  is  essential 
to  the  functioning  of  the  SOFC  and  cannot  be  done  away  with. 
While,  there  exits  a  possibility  that  the  decrease  in  the  partial 
pressure  of  the  reactants  due  to  poor  distribution  of  reactants 
along  the  length  of  the  fuel  cell  could  be  altered  with  a  differ¬ 
ent  design  of  the  flow  channel.  The  objective  of  this  paper  is  to 
explore  the  possibility  of  the  same. 

2.  Mathematical  model 

2.1.  Geometry  of  the  fuel  cell 

CFD  ACE  commercial  package  was  utilized  for  the  SOFC 
simulations.  The  CFD- ACE  package  discretizes  the  equations  of 
mass  momentum  and  energy  using  the  finite-volume  approach. 
SIMPLEC  (semi-implicit  method  for  pressure-linked  equations 
consistent)  [18],  a  variant  of  the  SIMPLE  algorithm  [19]  is  used 
to  solve  the  discretized  in  the  three-dimensional  space.  The 
porous  media  settings  in  the  package  utilizes  isotropic  linear 
resistance  model  (Darcy)  with  porosity  and  permeability.  The 
package  allows  for  the  electrochemical  (Butler- Volmer  kinet¬ 
ics)  heterogeneous  reactions  within  porous  media.  The  current 
continuity  equations  for  both  the  pore  phase  as  well  as  the  solid 
phase  of  porous  media  is  solved  in  the  package  which  facili¬ 
tates  in  the  accurate  determination  of  the  electrode  overpotential, 
necessary  to  treat  electrochemical  reactions  via  Butler- Volmer 
kinetics.  A  detailed  explanation  of  the  various  sub-models  of 
CFD  ACE  commercial  package  used  in  the  modeling  of  SOFC 
can  be  found  in  [20].  The  dimensions  of  the  co-flow  geometry 
SOFC  used  in  the  present  calculations  are  as  shown  in  Fig.  2.  The 
side  view  shows  the  arrangement  of  anode,  cathode,  electrolyte 
and  the  contacts. 
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Fig.  2.  Schematic  of  the  SOFC  used  in  the  calculations  along  with  its  various  parts  and  dimensions  (not  to  scale). 
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Table  1 


Porous  media  settings 


Volume  name 

£ 

K 

k 

RXN 

s/v 

Pore 

Diffusivity 

o 

Anode 

0.40 

le-12 

6.23 

Anode 

100 

le-6 

Bruggman(2) 

100000 

Cathode 

0.50 

le-12 

9.6 

Cathode 

too 

le-6 

Bruggman(2) 

7700 

Electrolyte 

0.01 

le-18 

2.7 

- 

- 

le-6 

Bruggman(2) 

le-20 

e:  porosity  (fluid  volume/total  volume);  k:  permeability  (m2)  (total  volume/surface  area)2;  k:  thermal  conductivity  (W  m  1  K  1 );  S/V:  surface  to  volume  ratio  (m  ); 
o:  electrical  conductivity  ( 51  ni  '  ). 


2.2.  Materials  and  property  settings 

The  anode  (Ni/ZrCh  ceramic  metal  composite  (cermet)), 
cathode  (doped  lanthanum  manganite  (LaMnOa))  and  elec¬ 
trolyte  (yttria-stabilized  Zr02  (8mol%  Y2O3))  are  modeled  as 
porous  media.  The  details  of  the  porous  media  settings  are  as 
given  in  Table  1 .  The  transport  properties  and  conductance  used 
in  the  current  calculations  are  as  given  in  Table  2.  Representa¬ 
tive  or  nominal  values  have  been  used  here  to  demonstrate  the 
viability  of  the  idea. 

2.3.  Chemical  reactions 

The  reactions  in  the  anode  and  the  cathode  are  modeled  as 
surface  reactions  and  are  assumed  to  take  place  inside  the  pores 
according  to  the  following  reaction.  The  oxygen  ions  are  mod¬ 
eled  as  “bulk  species”.  The  SOFC  model  used  here  does  not 
account  for  internal  reforming. 

anode  :  Hb  +  CO  +  202-  — >■  H2O  +  CO2  +  4e- 
cathode  :  O2  +  4e-  — >■  202- 

The  nominal  values  for  the  reference  current  Jo  and  the  Tafel 
constants  for  the  two  reactions  are  as  follows: 

Jq  —  1.E  +  14,  aa  =  0.7,  ac  =  0.7  for  anode  reaction 

Jo  —  IE  +  10,  aa  =  0.7,  ac  =  0.7  for  cathode  reaction 

2.4.  Boundary  conditions 

The  inlet  mass  flow  rate  into  the  cathode  and  the  anode  flow 
channels  were  2.12  x  10-6  and  3.415  x  10-7kgs-1,  respec¬ 
tively.  The  inlet  temperature  and  pressure  into  both  the  flow 
channels  were  1273  K  and  1  atm,  respectively.  At  both  the  cath¬ 


ode  and  anode  outlet,  the  fixed  pressure  boundary  condition  was 
used  and  the  pressure  was  1  atm. 

2.5.  Grid  structure  and  grid  details 

The  number  of  grid  points  utilized  in  the  present  study  is 
around  56,000.  The  typical  number  of  grid  points  in  the  x,  y 
and  z  direction  are  80,  70  and  10,  respectively.  The  results  were 
verified  to  be  grid  size  independent.  The  calculations  are  carried 
out  till  the  residuals  associated  with  all  variables  decrease  by  at 
least  4  orders  of  magnitude.  The  time  taken  for  convergence  is 
around  4  h  on  an  Intel  4  PC  with  0.5  GB  RAM  and  2.4  GHz. 

3.  Results  and  discussions 

The  modeling  of  the  flow  through  a  single  SOFC  was  car¬ 
ried  out  at  atmospheric  pressure  with  the  CFD-ACE  commercial 
package.  The  composition  of  the  fuel  used  was  0.096  H2  +  0.428 
H2O  +  O.26  CO +  0.216  CO2  and  air  composition  was  taken  as 
0.79  N2  +0.21  O2.  The  composition  of  the  fuel  is  a  typical  one 
obtained  from  an  externally  reformed  hydrocarbon  fuel.  In  order 
to  verify  the  computational  model,  the  variation  of  voltage  and 
power  density  with  the  current  density  (refer  Fig.  3)  was  plotted. 
The  results  show  a  qualitative  agreement  (similar  trend)  with 
those  reported  in  the  experimental  study  by  Jiang  and  Virkar 
[21],01gaetal.  [22]  and  Sasaki  et  al.  [23].  The  maximum  power 
from  the  SOFC  is  obtained  at  around  0.73  V.  The  calorific  value 
of  both  hydrogen  and  carbon  monoxide  is  nearly  the  same  and 
is  around  285  kJ  mole-1  (see  [1]).  Using  the  relationship  pre¬ 
sented  in  Larminie  and  Dicks  [16],  the  maximum  possible  EMF 
that  can  be  generated  from  the  SOFC  is  estimated  as  1 .48  V.  The 
efficiency  of  this  SOFC  at  maximum  power  condition  (refer 
Larminie  and  Dicks  [16])  is  around  50%. 

Thus  having  verified  the  computational  model,  a  thin  wall  as 
shown  in  Fig.  4  was  introduced  in  both  the  cathode  and  anode 


Table  2 

Transport  property  and  conductance  settings 


Volume  Name 

p 

P 

cr(  1  £2  1  m  *) 

CP 

k 

r 

Anode 

IGL 

MixKin 

10 

JANNAF 

MixKin 

SCH(0.7) 

Cathode 

IGL 

MixKin 

10 

JANNAF 

MixKin 

SCH(0.7) 

Electrolyte 

IGL 

MixKin 

10 

JANNAF 

MixKin 

SCH(0.7) 

Anode  channel 

IGL 

MixKin 

le-20 

JANNAF 

MixKin 

SCH(0.7) 

Cathode  channel 

IGL 

MixKin 

le-20 

JANNAF 

MixKin 

SCH(0.7) 

p:  density  (kgm  3);  IGL:  ideal  gas  law;  p:  viscosity  (kgm  1  s  *);  MixKin:  kinetic  theory;  cr:  electrical  conductivity  (£2  1  m  *);  Cp :  specific  heat  (Jkg  1  K  *)  by 
JANNAF  curve  fits;  k\  thermal  conductivity  (W m-1  K-1)  and  T"1:  mass  diffusivity  (kgm-1  s_1)  by  Schmidt  number. 
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Fig.  3 .  The  polarization  characteristics  and  the  power  density  of  a  plain  geometry 
SOFC. 

flow  channel.  The  purpose  of  introducing  such  a  thin  wall  as 
already  been  explained  in  the  “Introduction”  section  is  to  have  an 
equitable  distribution  of  reactants  along  the  length  of  the  SOFC. 
It  is  seen  in  Fig.  4  that  the  thin  wall  runs  up  to  half  the  length 
of  the  fuel  cell  and  that  it  divides  the  inlet  sections  (anode  and 
cathode  flow  channel)  into  two  halves.  Thus,  fresh  reactants  are 
also  made  available  to  the  second  half  (Section  II,  refer  Fig.  4) 
leading  to  an  even  distribution  of  fresh  reactants  over  the  two 
half  lengths  (Sections  I  and  II,  refer  Fig.  4)  of  the  fuel  cell. 

For  the  new  geometry  of  the  SOFC,  the  calculations  were  car¬ 
ried  out  to  obtain  the  polarization  characteristics.  The  inlet  mass 
flow  rates  and  all  other  conditions  were  kept  identical  for  both 
the  plain  geometry  and  the  thin-walled  geometry.  The  results 
of  these  calculations  for  both  the  plain  geometry  and  the  thin- 
walled  geometry  are  as  illustrated  in  Fig.  5.  This  figure  shows  the 
variation  of  voltage  and  power  density  with  the  current  density.  It 
is  seen  from  Fig.  5  that  at  all  voltages  the  current  density  obtained 
from  the  thin-walled  geometry  is  higher  than  those  obtained  with 
the  plain  geometry.  The  improvement  in  the  output  power  with 
the  thin-walled  geometry  over  the  plain  geometry  can  also  be 
observed  from  Fig.  5.  The  peak  power  is  obtained  with  the  thin- 
walled  geometry  at  around  0.88  V,  whilst  the  same  with  a  plain 
geometry  is  obtained  at  0.73  V.  Thus,  the  efficiency  of  the  SOFC 
at  peak  power  for  the  thin-walled  geometry,  proceeding  along 


Current  Density.  A/cm2 

Fig.  5.  The  polarization  characteristics  and  the  power  density  of  both  plain  and 
thin-walled  geometry  SOFC. 

lines  similar  to  those  outlined  in  conjunction  with  Fig.  3  is  60%. 
The  peak  power  density  obtained  with  the  thin-walled  geometry 
is  1.18Wcm-2  and  the  same  obtained  with  a  plain  geometry 
is  0.88  W  cm~2.  The  peak  power  density  is  around  35%  higher 
with  the  thin-walled  geometry  compared  to  the  plain  geometry. 
Thus,  with  the  introduction  of  the  thin  wall  in  the  flow  chan¬ 
nel  as  shown  in  Fig.  4,  both  the  efficiency  and  the  maximum 
power  density  of  the  SOFC  were  enhanced.  Here,  it  must  be 
emphasized  that  the  relative  improvement  in  performance  of  the 
thin-walled  geometry  over  the  plain  one  is  more  important  than 
the  quantitative  values  of  power  density  obtained  with  the  two 
geometries.  This  method  will  work  when  one  is  constrained  to 
use  air  and  combination  of  different  gases  as  fuel  as  is  currently 
used.  If  hydrogen  is  used  as  fuel  and  oxygen  is  used  as  oxidizer, 
the  introduction  of  thin-walled  geometry  will  not  enhance  the 
performance  significantly  (as  reported  here)  due  to  the  absence 
of  diluants.  But  in  many  practical  applications,  due  to  cost  com¬ 
pulsions  and  other  considerations  air,  instead  of  oxygen  and  a 
combination  of  fuels  instead  of  single  fuel  is  used.  The  intro¬ 
duction  of  thin-walled  geometry  in  such  cases  could  enhance  the 
performance. 

To  understand  the  processes  that  led  to  the  better  performance 
of  the  thin- walled  geometry  compared  to  the  plain  geometry,  the 
current  density  variation  along  the  axis  (x)  at  y  =  0  and  z  =  0  for 
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Fig.  4.  Schematic  of  the  thin  walled  geometry  SOFC  used  in  the  calculations  along  with  the  mass  flow  rates  at  various  locations.  The  dimensions  are  as  indicated 
(not  to  scale). 
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Fig.  6.  The  current  density  variation  along  the  axial  length  of  both  the  plain  and 
thin-walled  geometry  SOFC  at  y  =  0  and  z  =  0  at  a  voltage  of  0.73  V. 

the  two  geometries  were  plotted  (refer  Fig.  6).  The  current  yield 
for  both  geometries  from  ‘Section  I’  (refer  Fig.  4)  is  almost 
identical  as  the  area  under  the  two  curves  are  almost  the  same. 
The  current  yield  from  the  second  half  or  ‘Section  II”  (refer 
Fig.  4)  of  the  SOFC  with  the  thin-walled  geometry  is  higher 
compared  to  the  plain  geometry  SOFC.  Similar  trends  were 
observed  at  other  (y  and  z)  locations  too.  The  increase  in  the 
current  density  in  Section  II”  (refer  Fig.  4)  of  the  SOFC  with 
the  thin-walled  geometry  can  be  attributed  to  the  availability  of 
fresh  reactants  (explanation  to  follow  later  in  the  paper).  Inci¬ 
dentally,  if  the  plain  geometry  were  also  to  be  divided  into  two 
halves  along  its  length,  then  the  power  density  obtained  from 
it’s  first  (nearer  to  inlet)  and  second  (nearer  to  outlet)  halves  are 
0.75  and  0.4425  W  cm-2,  respectively.  The  corresponding  val¬ 
ues  for  the  thin-walled  geometry  are  0.8 125  and  0.6375  W  cm-2, 
respectively.  The  performance  of  the  thin-walled  geometry  is 
slightly  higher  than  the  plain  geometry  in  the  first  half,  due  to 
the  increased  mass  flux  made  possible  in  the  thin- walled  geome¬ 
try.  The  performance  of  the  thin- walled  geometry  is  significantly 
higher  than  the  plain  geometry  in  the  second  half.  This  indicates 
that  the  overall  enhancement  in  performance  observed  with  the 
thin-walled  geometry  is  probably  due  to  the  better  distribution 
of  reactants  to  Section  II. 

The  value  of  current  density  at  x  =  0  for  both  geometries  are 
different  from  each  other  as  can  be  observed  from  Fig.  6.  This 
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Fig.  7.  The  pressure  distribution  along  the  length  of  the  SOFC  close  to  anode 
and  cathode  for  both  plain  and  thin-walled  geometry  at  a  voltage  of  0.73  V. 


is  a  rather  strange  behavior  as  the  initial  composition  of  the  fuel 
and  air  is  identical  for  both  geometries.  To  understand  the  reason 
behind  this  the  pressure  plots  close  to  the  anode  (y=  —0.525  mm) 
and  cathode  (y  =  0.525  mm)  were  plotted  (refer  Fig.  7).  It  is  evi¬ 
dent  from  Fig.  7  that  the  pressure  increases  with  the  introduction 
of  the  thin  wall  in  both  anode  and  cathode.  From  Darcy’s  law, 
an  increase  of  pressure  near  the  anode  and  cathode  (both  porous 
media)  means  increased  mass  diffusing  into  anode  and  cathode. 
This  increased  mass  flow  leads  to  greater  reactions  and  hence 
higher  current  density  (atx  =  0)  in  the  thin- walled  geometry. 

In  order  to  understand  the  cause  for  the  improvement  in  per¬ 
formance  of  the  thin-walled  geometry  over  the  plain  geometry, 
the  mass  fraction  profiles  of  a  representative  species  CO  (the 
largest  reactive  fuel  fraction,  CO  26%  and  FF  around  10%)  were 
plotted  for  both  the  plain  geometry  and  the  thin-walled  geometry. 
They  are  as  shown  in  Fig.  8.  It  is  seen  that  with  the  introduction 
of  the  thin  wall  in  the  flow  channel,  the  distribution  of  reactants 
along  the  length  of  the  flow  channel  improves  significantly.  On 
the  basis  of  Fig.  8,  it  can  be  inferred  that  the  improvement  in  the 
performance  of  the  thin-walled  geometry  over  the  plain  ones  is 
due  to  the  more  even  distribution  of  reactants  along  the  length 
of  the  SOFC. 


Plain  Geometry 


>  Y 


Fig. 


The  mass  fraction  profiles  of  CO  for  both  the  plain  and  the  thin-walled  geometry  (Vertical  lines  are  an  artifact  of  the  griding  procedure  only). 
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Table  3 

Current  density  yields  for  the  thin  walled  and  plain  geometries  at  a  voltage  of  0.73  V 


383 


Factor  by  which  the  basic 
flow  has  been  increased 

Current  density  for  the  plain 
geometry,  (A  cm-2) 

Current  density  for  the 
thin- walled  geometry,  (A  cm-2) 

Improvement  of  current  density 
with  thin-walled  geometry  (%) 

0.5 

0.805 

0.8625 

7 

1 

1.1925 

1.4537 

21.9 

2 

1.5025 

1.8538 

23.4 

4 

1.5962 

1.9850 

24.4 

6 

1.6187 

2.0200 

24.8 

8 

1.6275 

2.0375 

25.2 

10 

1.6362 

2.0488 

25.2 

The  enhancement  in  the  performance  of  the  SOFC  with  the 
thin-walled  geometry  at  a  particular  inlet  flow  rate  of  reactants 
has  been  observed.  But  is  this  enhancement  in  performance 
present  at  different  inlet  mass  flow  rates  as  well?  To  address 
the  above  question,  calculations  on  the  thin-walled  geometry 
and  the  plain  geometry  were  carried  out  at  different  inlet  mass 
flow  rates.  The  results  of  these  calculations  are  as  tabulated  in 
Table  3  for  a  fixed  voltage  of  0.73  V.  The  cathode  and  the  anode 
mass  flow  rates  of  2. 12  x  10_6and3.415  x  1CP7  kgs-1,  respec¬ 
tively  were  used  as  the  basic  unit.  The  first  column  in  Table  3 
gives  the  factor  by  which  the  basic  cathode  and  anode  inlet  mass 
flow  rates  have  been  simultaneously  altered.  The  last  column  in 
Table  3  shows  that  the  performance  enhancement  obtained  with 
thin- walled  geometry  over  the  plain  geometry.  The  performance 
enhancement  with  the  thin-walled  geometry  varies  from  7-25% 
over  the  plain  geometry. 

At  inlet  mass  flow  rates  lower  than  the  basic  unit  (cathode  and 
anode  mass  flow  rates  of  1.06  x  10-6  and  1.707  x  10~7  kg  s-1, 
respectively)  low  levels  of  improvement  (7%)  are  observed.  At 
these  low  inlet  mass  flow  rates,  the  quantity  of  reactants  are 
lower  and  they  get  greater  residence  time  to  react  as  the  flow 
velocities  are  halved.  Hence,  a  small  enhancement  (7%)  in  per¬ 
formance  is  obtained  with  the  introduction  of  thin  wall.  As  the 
inlet  mass  flow  rates  are  increased  to  the  basic  unit  and  above  it, 
the  current  densities  obtained  from  the  system  are  also  increased 
(refer  column  2  and  3  of  Table  3).  This  is  similar  to  the  results 
reported  in  literature  by  Leng  et  al.  [24]  and  Larrain  et  al.  [25]. 
At  higher  inlet  mass  flow  rates,  higher  current  densities  (refer 
column  2  and  3  of  Table  3)  and  in  turn  higher  power  densities  are 
obtained.  Thus,  higher  inlet  mass  flow  rates  are  to  be  preferred, 
if  enhancing  the  power  density  is  the  objective.  At  high  inlet 
mass  flow  rates,  the  reactant  quantity  is  larger  and  the  accom¬ 
panied  residence  time  to  react  is  lower.  These  conditions  are 
quite  stringent  and  call  for  better  flow  management  of  reactants. 
Hence,  with  the  introduction  of  the  thin  wall  at  high  inlet  mass 
flow  rates,  the  improvement  in  performance  of  SOFC  obtained 
is  significant  (22-25%). 

4.  Conclusions 

Three-dimensional  simulation  of  a  single  SOFC  was  carried 
out  for  two  geometries  (plain  and  the  thin-walled  geometry) 
using  CFD-ACE  package.  The  thin-walled  geometry  had  a  thin 
wall  that  split  the  inlet  section  of  both  the  anode  and  cathode  flow 
channel  into  two  parts  and  ran  up  to  half  the  length  of  the  SOFC. 


All  conditions  remaining  identical  for  the  two  geometries: 

(1)  The  SOFC  performance  was  better  with  the  thin-walled 
geometry  as  compared  to  the  plain  one  due  to  the  better 
distribution  of  reactants  along  the  length  in  the  thin-walled 
geometry. 

(2)  The  thin-walled  geometry  not  only  yielded  a  higher  max¬ 
imum  power  density  of  1.18  W  cm-2,  but  it  also  led  to  an 
improvement  in  the  efficiency  at  which  the  maximum  power 
density  was  obtained  (around  60%).  For  the  plain  geometry 
the  maximum  power  density  of  0.88  W  cm-2  was  obtained 
at  an  efficiency  of  50%.  The  peak  power  density  is  around 
35%  higher  with  the  thin-walled  geometry  as  compared  to 
the  plain  geometry. 

(3)  The  improvement  in  the  performance  of  the  thin-walled 
geometry  over  the  plain  geometry  was  small  (7%)  at  low 
inlet  mass  flow  rates  and  increased  (~25%)  with  the  increase 
in  inlet  mass  flow  rates. 

This  research  opens  up  new  vistas  to  explore  the  possibility  of 
increasing  the  power  density  of  a  SOFC.  The  results  presented 
here  can  perhaps  be  viewed  as  a  logical  extension  of  the  desire 
to  have  a  uniform  flow  in  the  fuel  cell  stacks  in  order  to  improve 
its  effectiveness.  Although,  the  results  presented  here  are  for  a 
SOFC,  they  are  more  generic  in  nature  and  can  be  applied  to 
any  other  kind  of  fuel  cell.  The  ideas  presented  here  need  to  be 
experimentally  authenticated  before  proceeding  further. 
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